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Abstract 

Calculations of the chiral condensate (\x) on the lattice using staggered 
fermions and the Lanczos algorithm are presented. Four gauge fields are 
considered: the quenched non-Abelian field, an Abelian projected field, and 
monopole and photon fields further decomposed from the Abelian field. 
Abelian projection is performed in maximal Abelian gauge and in Polyakov 
gauge. The results show that monopoles in maximal Abelian gauge largely 
reproduce the chiral condensate values of the full non-Abelian theory, in both 
SU(2) and SU(3) color. 
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I. INTRODUCTION 



Since the Abelian monopole mechanism for confinement in QCD was first proposed [[!],§ , 
there have been extensive studies in the pure gauge sector of lattice theories The effect 



of Abelian projection in the quark sector was studied in where chiral symmetry breaking 
and meson correlators were analyzed. Further studies of the role of Abelian monopoles in 



chiral symmetry breaking at finite temperature [11|, and in hadron spectroscopy [12|, have 
been made. These works have provided evidence in support of Abelian monopoles as the 
mechanism of chiral symmetry breaking in QCD. 

In the present work, a systematic study of the role of Abelian projection and monopoles 
in chiral symmetry breaking is carried out in the confined phase at zero temperature. The 
chiral condensate is computed using four sets of gauge fields: the quenched non-Abelian 
field, an Abelian projected field, and monopole and photon fields further decomposed from 
the Abelian projected field. Abelian projection is performed in maximal Abelian gauge and 
in Polyakov gauge. We employ a different technique to compute the chiral condensate than 
was used in Ref. |[L0| , where an extrapolation to the chiral limit was made from nonzero 
quark mass: here we use the spectral representation of the chiral condensate along with the 
Lanczos algorithm which allows calculations to be done directly at zero mass. Most 
of the results presented here are for SU(2) color. Some results of a first study of chiral 
symmetry breaking in the Abelian projection of SU(3) gauge theory are also given. 

The calculational method is described in Sec. || and results are presented in Sec. fU . 
Using the Lanczos method at zero quark mass we find reasonable quantitative agreement 
between the chiral condensate calculated with the monopole part of the Abelian projected 
field and the quenched non-Abelian field for the SU(2) theory. Our present SU(3) simulation 
is not yet good enough to draw the same quantitative conclusion, but qualitatively, the same 
pattern emerges as seen in the SU(2) calculation. This provides some evidence that the 
Abelian dominance idea can be extended to the SU(3) theory. 



Some of our results were presented in summary form in Ref. [14 



II. METHOD 
A. Chiral Condensate on the lattice 

Chiral symmetry breaking is studied using staggered fermions, with the action 

S f = o H ^( x ) \x( x ) U i*( x ) x{% + A) ~ X& + fr) U t( x ) X{x)\ + m x(x)x(x), 

= XM({U}) X , (1) 

where x-, x are single- component fermion fields, rj^x) is the staggered fermion phase ||15|| , m 
is the mass in lattice units and the U's are gauge field links. 

The chiral symmetry order parameter is calculated from the inverse of the fermion matrix 
M of Eq. (D 

{xx(m,V)) = ±(TrM-\{U})), (2) 
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where V is the lattice volume and the angle brackets denote the gauge field configuration 
average. This can be rewritten in a spectral representation [16 



— 1 N 1 —1 2m 

fe(m ,^) ) = _|;_ = (3) 

where iX n 's are the eigenvalues of the zero-mass fermion matrix (Dirac operator). The 
eigenvalues appear in complex conjugate pairs so only the positive half of the spectrum 
needs to be considered as in the second part of Eq. (|3]). 

To correctly probe the physics of spontaneous chiral symmetry breaking, one should 
attempt to work in the limit of zero quark mass and infinite volume. The chiral limit m — ► 
should be taken after V — > oo 

(XX) = - lim lim I £ -^— 2 = - lira H dX^{ = -vrp(O), (4) 
m->o v^oo V A > A„ + m z m->o Jo A^ + m z 

where the spectral density function p(X) = ^dn/dX is normalized to / °° dXp(X) = N c , 
the number of colors. Eq. relates chiral symmetry breaking to the small modes in the 
eigenvalue spectrum. So the task is reduced to finding the small eigenvalues of the fermion 
matrix at zero mass, rather than the entire spectrum. This can be done using the well- 



established Lanczos algorithm [13 



B. Abelian Projection on the Lattice 

The lattice formulation of Abelian projection was developed in The idea is to fix the 
gauge of a SU(N) theory so that a residual [/(l)^ -1 gauge symmetry remains. The Abelian 
degrees of freedom are extracted by a subsequent projection U(x, p) = c(x, p) u(x, p), where 
u(x, p) is the (diagonal) Abelian-projected field and c(x, p) the non-diagonal matter field. 
In general the gauge condition can be realized by diagonalizing some adjoint operator 1Z: 

G(x)TZ(x)G~ 1 (x) = diagonal. (5) 

Several gauge conditions have been studied and it has been found that the so-called 
maximal Abelian gauge f| most readily captures the long-distance features of confinement 
(the relevance of other gauges to long-distance physics has been considered in 0). In SU(2) 
this gauge is realized by maximizing the quantity 

R = J2te[°3U(x,p)(T 3 W(x,p)] (6) 

where U(x, p) = G(x)U(x, p)G~ 1 (x+p). We maximize R by an iterative procedure. We find 
a gauge transformation G(x) that maximizes R locally, at given site x, by keeping G(x + p) 
at neighboring sites fixed. Using the fact a^G = GV3 (with the exception of diagonal gauge 
transformations) leads to an analytical solution for the local maximization of R 

G 2 (x)=F-\x), (7) 
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where 

4 

F ( x ) = E [u(x,Li)a 3 U Jf (x,n)a 3 + U j (x - n,n)a 3 U(x - n,n)a 3 . (8) 



The global maximum of R is approached iteratively by repeatedly sweeping through the 
lattice, until G(x) is sufficiently close to the identity: 

max jl - iTrG(x) j < S < 1, 5 ~ 1(T 6 . (9) 

Maximal Abelian gauge in SU(3) is implemented by maximizing 

R = E [\Uu(x^)\ 2 + (^(x,^)! 2 + |t/ 33 (^/i)l 2 ] • (io) 

The maximization is done by going through the three SU(2) subgroups of SU(3). The gauge 
transformation within the first subgroup can be parameterized as 



go 92 + igi o > 
Gi{x) = | -g 2 + m go 
1/ 



(11) 



with the constraint + g\ + g\ = 1. R is maximized locally at x, again by holding G± at 
neighboring sites fixed. Under this local transformation, R can be expressed as 

R = E [°i9o + c 29i + C3#2 + dgogi + c 5 g g 2 ] (12) 

x 

where the sum over /i is implicit in the coefficients q. The solution is obtained nu- 
merically, and this is facilitated by introducing spherical coordinates on the unit sphere: 
#2 = sin 9 cos 0, gi = sin 9 sin 0, g = cos 9. This procedure is repeated at a given site for 
the second and the third subgroups. The global maximum of R is approached by repeatedly 
sweeping through the entire lattice until a stopping criterion similar to Eq. (|9|) is satisfied. 

Once the gauge fixing is done, the Abelian-projected links are extracted. In SU(2) gauge 
theory, where the links can be parameterized by U(x, /i) = r + ia ■ f, the Abelian projected 
links u are given by 

u = diag(e^, e"^), <p = tan" 1 ^/^). (13) 
In SU(3), the Abelian configurations are extracted according to 

u = diag(tti,tt 2 ,u 3 ), «fc = exp (i axgU kk - ^<p), (14) 

where 

V 9 = (^2 aT SU k kj mod 2tt C (-7r,7r]. (15) 

The three phase factors are constrained by U\U2U 3 = 1 so that only two of them are inde- 
pendent. 
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C. Monopole Decomposition 



It is well-known that there exist monopoles in a compact U(l) field. The Abelian- 
projected link u(x,fi) = exp(i0 M (a:)) can be resolved into a component due to monopoles 
by considering the Abelian field strength (fi^x), defined from plaquette phases in the usual 
way 

</vO) = d^(p u (x) - d u (p^x), (16) 

where d^f(x) — f(x + /t) — f(x). The flux due to an integer- valued monopole string m^ v is 
identified from the field strength according to [[IS 



;>vk x ) = ^»( x ) + 2vrm^(x), (17) 



where <\>' G (—71,71]. The vector potential 0™ on (a;) generated by the monopoles is therefore 
given (in Lorentz gauge) by 

C on (x) = -2nJ2D(x-y)d^m, u (y), (18) 

y 

where the lattice photon propagator satisfies —d , d IJl D(x) = S Xt Q, with & f(x) = f(x) — f(x — 
fi) [[H]]. The "photon" field 0^ is identified with the difference M (x) 



b™ OIi (x) 



III. NUMERICAL RESULTS 



The Lanczos algorithm |13| is carried out in double precision to find the eigenvalues 



of the fermion matrix. To test our implementation of the algorithm, the full spectra of 
an 8 lattice and a 12 lattice at (3=2.3 in SU(2) were obtained for one configuration. In 
SU(2) one knows a priori that every eigenvalue is doubly degenerate so that the number 
of single positive eigenvalues must be V/2. Furthermore, they satisfy the closure relation 
J2n=i = V. We found exactly 2048 single positive eigenvalues for the 8 4 lattice and 
10368 for the 12 4 lattice, and the closure relations were satisfied to a few parts in 10 8 . This 
indicates that we can determine the spectrum accurately. 

A. Results for SU(2) 

A heat-bath Monte Carlo algorithm was used to generate quenched gauge field configu- 
rations using the standard Wilson plaquette action, with periodic boundary conditions on 
a 14 4 lattice at /3=2.2, 2.3, 2.4 and 2.5. Anti-periodic boundary conditions were used for 
fermions in all directions. Other boundary conditions were also tried, and our results showed 
little change. Gauge fixing was done with the help of overrelaxation [HI which reduced the 



number of iterations by a factor of 3 to 5. About 500 iterations with overrelaxation were 
required for a stopping criterion 5 ~ 10 -6 . 

In Fig. we show the raw data for the spectral density function p(X) obtained from 70 
configurations at each j3 value. The results clearly show a non-vanishing signal of p(A) at all 



5 



(3 values, although finite volume effects manifested as a sudden depletion near A = pi] 



begin to set in at j3=2A and 2.5. To extract a value at A = 0, we fit the distributions by a 
straight line p(A) = p(0) +p'(0)A in an interval [A mi „, A max ]. The interval is chosen so that it 
excludes those eigenvalues near A = that are strongly influenced by finite volume effects, 
and those that cause p(A) to depart from linear behavior. It is expected that X m in ~ £, 3 /V 
where £ is some length scale governed by the gluon dynamics [^TJ. The value A max is chosen 
to be as large as possible while preserving both the stability and the quality of the fit. Under 
these conditions, we find that the fit is quite stable over a relatively wide interval. Fig. |2| 
shows the results of such a fit for the eigenvalue interval [0.015,0.05]. The fitted values for 
p(0) and p'(0) and the extracted chiral condensate (xx) are given in Table | in lattice units. 
The errors quoted are statistical and are obtained using the jackknife method. Since the 
eigenvalue spectrum is doubly degenerate, the quoted numbers for p(A) are only half of the 
true value, so that the extracted (xx) = — 27rp(0). The value of (xx) i n the full SU(2) theory 
and in maximal Abelian gauge projection are consistent with those obtained in Ref. |TU 



where an extrapolation to the chiral limit was made from nonvanishing quark mass. The 
interesting feature here is that the monopole configuration average has a condensate that is 
even closer to the full theory than the Abelian projected fields. On the other hand with the 
photon configurations either no or very few small eigenvalues are found. Hence we conclude 
there is no chiral symmetry breaking from these configurations. 

For purposes of comparison, we also performed Abelian projection at (3 — 2.5 using a 
different gauge-fixing condition: the Polyakov gauge, in which the Polyakov loop is diago- 
nalized according to Eq. (^|). The result is shown in Fig. |3|. We see that the Abelian and 
the monopole field spectral density functions are almost an order of magnitude larger than 
those of the full theory. In Ref. a similar result for (xx) * s found using the field-strength 



gauge. 



B. Results for SU(3) 

The gauge field configurations were generated using the Cabibbo-Marinari J20| pseudo 



heat-bath method on a 8 4 lattice at (3=5.5, 5.6, 5.7, 5.8 and 10 4 lattice at (3 = 5.9. Configu- 
rations are selected after 4000 thermalization sweeps from a cold start, and every 500 sweeps 
thereafter. Fig. f| shows the spectral density function obtained from 44 SU(3) gauge field 
configurations for (3=5.5, 66 configurations for (3=5.6, 150 configurations for (3=5.7, and 37 
configurations for (3=5.8. 

We performed Abelian projections on the 8 4 lattice at /3 = 5.7 and 10 4 lattice ai (3 = 5.9. 
Gauge fixing in SU(3) is time-consuming. On average about 500 iterations with overrelax- 
ation are required for 5 = 1 — l/9Tr[Gi(x) + G^O^) + Gs(x)] to converge to 10~ 5 . After 
configuration averaging, the three phase factors in Eq. ([i~4[) give equal contributions to the 
chiral condensate. So the total is obtained by first calculating with one phase and then 
multiplying by 3. This saves a factor of 3 in computer time. Fig. |5] and Fig. |6] show the data 
obtained for p(A). For the 8 4 lattice at (3 = 5.7, 150 configurations were used, and 100 con- 
figurations for 10 4 at /3 = 5.9. We see that a similar pattern emerges in SU(3) as in SU(2): 
for small eigenvalues, the Abelian and the monopole contributions give spectral densities 
that are close to those of the full theory. It was also confirmed that photon configurations 
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give negligible contributions. 



IV. CONCLUSION 

We have calculated chiral symmetry breaking on the lattice in the quenched approxima- 



tion and Abelian projection. Unlike the previous study HlOU , the Lanczos method is used to 
calculate eigenvalues of the fermion matrix and hence the chiral condensate directly at zero 
quark mass. Furthermore, in this work, a decomposition of the Abelian projected field into 
monopole and photon pieces was made. For SU(2) gauge theory it was found the monopole 
part of the Abelian field, projected in maximal Abelian gauge, yields chiral condensate val- 
ues which are quite close to those obtained with the full non- Abelian fields. In contrast, the 
photon piece of the Abelian projected field gives no condensate. 

A calculation was also done using Abelian fields projected in the Polyakov gauge. This 
yielded an eigenvalue density and hence a chiral condensate about an order of magnitude 
larger than the non- Abelian calculation. This is consistent with what was previously found 



in field strength gauge in Ref. |TU 



Some calculations were also done in an SU(3) theory. Qualitatively, a similar pattern is 
seen as in the SU(2) calculation. In the region of small eigenvalues (which is relevant for 
chiral symmetry breaking), the Abelian and monopole fields give spectral densities which 
are quite close to that of the full non- Abelian calculation. This provides a positive indication 
that the idea of Abelian dominance can be extended to the SU(3) theory. 
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FIG. 1. Spectral density functions calculated at four different (3 values for three types of 

gauge configurations: the non-Abelian SU(2) configurations, the Abelian-projected configurations 

in the Maximal Abelian gauge, and the monopole configurations further decomposed from the 

Abelian-projected configurations. 
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FIG. 2. Fitted spectral density functions in the interval [0.015,0.05] in SU(2). 
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FIG. 3. Spectral density functions calculated in the Polyakov gauge at (3=2.5 for three types 
of gauge configurations: non-Abelian SU(2)(»), Abelian(D) and monopole(V). 



10 





FIG. 4. Spectral density functions calculated at four different (3 values in SU(3). 
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TABLES 



TABLE I. The extracted values for 


(xx) 


on the 14 4 lattice from the straight line fit 


p(A) = p(0) + p'(0)X in the interval [0.015,0.05] 


At each value, 


the first row is for the full 


SU(2), second row Maximal Abelian, third 


row monopole. Numbers 


in parentheses are statistical 


errors in the final digit (s). 
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